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1.  INTRODUCTION 


This  report  describes  the  results  of  computer  calculatiorts  to  find  the  pressure  distributions  on  two  axisymmetric  bomb 
shapes  of  current  interest  to  Aerodynamic  Research  Group  at  W.R.E.  The  studies  were  initiated  to  augment  wind-tunnel 
tests  of  the  two  bodies,  and  were  particularly  concerned  with  the  effects  of  changes  to  the  base  geometry,  including  the 
influence  of  the  sting  on  wind-tunnel  models. 

The  calculations  were  limited  to  inviscid,  axisymmetric  flow  solutions,  and  because  the  emphasis  was  on  simplicity  and 
case  of  operation,  the  method  chosen  was  that  of  Landweber(ref.l)  which  has  recently  been  programmed  in  Fortran  by 
Albone(ref.2).  The  theoretical  formulation  assumes  incompressible  flow  and  the  results  given  here  are  for  zero  Mach 
number;  however  a Gothert-type  correction  has  been  included  to  enable  approximate  solutions  for  subsonic  compressible 
flow  to  be  obtained. 


2.  DESCRIPTION  AND  USE  OF  THE  PROGRAMS 


2.1  The  integral  equation 

The  method  is  equivalent  to  a representation  of  the  body  by  a ring  vortex  distribution  r(s)  over  its  surface. 
With  the  coordinate  system  shown  in  figure  1 , the  velocity  induced  at  station  t on  the  body  axis  by  the  vorticity 
on  the  surface  element  ds  is  (see,  for  example,  reference  3); 


(x)r(s)ds 

U(x,t)=-^ . (1) 

2[y^(x)  + (x-ty  i 


A consequence  of  representing  the  body  surface  by  a vortex  sheet  is  to  replace  the  body  interior  by  fluid  at  rest. 
Thus  the  vortex  distribution  must  induce  an  interior  velocity  -U»  to  cancel  the  incident  free  stream.  In  particular, 
the  velocity  induced  along  the  axis  must  be  -U,» . 

It  can  be  shown  that  the  velocity  jump  across  a vortex  sheet  is  minus  the  local  vortex  intensity.  The  condition 
of  zero  total  internal  velocity  then  gives 


u(x)  = - P(s), 


(2) 


where  u(x)  is  the  longitudinal  surface  velocity  at  station  x.  Substituting  (2)  into  (1)  and  integrating  along  the 
body  surface  gives  the  total  velocity  on  the  body  axis,  and  equating  to  -U<»  gives  the  integral  equation; 

r ^ u(x)  y^  (x)  ds 

2 — = U.  (3) 

\ 2[y^(x)  + (x-t)M  ^ 

to  be  satisfied  at  all  points  t between  O and  P.  The  solution  of  (3)  gives  the  velocity  distribution  u(x),  and  thence 
the  pressure  distribution  Cp  (x). 

A rigorous  derivation  of  (3)  and  an  iterative  method  of  solution  have  been  given  in  reference  1 and  repeated  in 
reference  2,  which  also  describes  Fortran  programs  to  solve  for  both  closed  bodies  and  those  with  a parallel 
afterbody  extending  downstream  to  infinity.  These  programs  need  only  minor  alterations  to  be  run  on  the  W.R.E. 
computing  system.  For  intending  users  without  ready  access  to  reference  2,  listings  are  included  in  the  Appendix. 

The  fact  that  the  integral  equation  (3)  for  the  surface  velocity  is  satisfied  on  the  body  centreline,  means  that 
the  method  may  not  be  as  accurate  as  others  which  satisfy  a boundary  condition  on  the  body  surface,  such  as 
that  of  Hess  and  Smith(ref.4).  For  body  shapes  not  too  dissimilar  from  an  ellipsoid  of  revolution,  the  solution 
converges  to  any  desired  degree  of  accuracy,  but  for  more  complicated  shapes  the  solution  may  converge  to  a best 
solution  and  then  become  divergent.  In  such  cases  this  best  solution  is  usually  sufficiently  accurate  for  most 
practical  purposes.  In  any  case,  the  programs  are  simple  to  use  and  require  only  minimal  computing  times 
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coinparcd  with  the  more  accurate  methods,  which  determined  their  use  in  the  present  application.  Another 
property  of  the  method  of  solution  is  that  corners  or  sudden  changes  of  curvature  in  the  body  profile  tend  to  be 
“smoothed”,  which  is  not  a serious  limitation  since  the  boundary  layer  will  have  a similar  effect  in  the  real  flow. 

2.2  Input  and  output 

Input  to  each  program  consists  of  a subroutine  called  BODY  which  can  calculate  the  radius  and  surface  slope  of 
the  body  at  any  longitudinal  station,  and  several  parameters  which  are  declared  as  data  at  the  beginning  of  the 
program.  For  the  closed-body  program  these  are  as  follows; 

NMAX  - the  maximum  number  of  iterations  allowed  (e.g.,  200). 

EPS  - a quantity  which,  when  greater  than  the  maximum  difference  between  any  two  successive 

iterations,  causes  the  iteration  procedure  to  stop. 

XXO  - abscissa  of  body  nose. 

XXI  - abscissa  of  body  tail. 

MACH  - Mach  number,  used  with  a Gothert-type  correction  to  allow  for  compressibility  effects. 

For  the  infinite  - afterbody  program,  the  input  parameters  are  as  follows: 

NMAX,  EPS,  MACH  - as  for  the  closed-body  program. 

LENGTH  - length  of  forebody.  Note  that  in  this  program,  BODY  assumes  the  forebody  to  be 
between  x = 0 and  x = LENGTH,  so  that  for  complicated  bodies  the  same  basic 
subroutine  can  be  used  for  each  program. 

YAFT  - radius  of  afterbody. 

The  program  output  consists  of  a listing  of  the  Gaussian  abscissae  used  in  the  numerical  integration  for  each 
iteration  (40  points  along  the  body  axis  in  the  case  of  the  closed-body  program,  and  30  points  along  the  forebody 
axis  for  the  infinite-afterbody  program),  together  with  the  associated  values  of  body  radius  and  slope,  non- 
dimensional  surface  velocity  and  pressure  coefficient.  The  number  of  iterations,  and  maximum  difference  between 
the  two  final  iterations,  are  also  printed. 

2.3  Body  geometry 

Two  basic  configurations  of  interest  were  studied,  and  designated  Body-A  and  Body-B  respectively.  Two  further 
variations  of  the  Body-A  tail  geometry  are  referred  to  as  A(i)  and  A(ii).  The  half-body  profiles  of  the  various 
configurations  are  given  in  figures  2 to  6,  the  x - values  being  in  body  calibres. 

In  those  cases  where  the  configuration  has  a blunt  base,  for  example  when  modelling  the  free-flight  situation, 
some  extension  of  the  model  surface  is  necessary  to  simulate  the  surface  streartdine  separating  from  the  base. 

These  extensions  are  shown  on  the  figures  as  a broken  line.  Where  a sting  whose  diameter  is  equal  to  the  base 

diameter  d„  is  present  this  problem  does  not  arise  and  the  sting  is  assumed  to  extend  downstream  to  infinity. 

B 


3.  RESULTS  AND  DISCUSSION 


3.1  Standard  Body-A  configuration 

Figure  2 shows  the  variation  of  surface  pressure  coefficient  on  the  standard  Body-A  geometry  with  base  flow 
represented  by  (a)  a strong  flow  expansion  behind  the  base;  (b)  a large  sting  of  diameter  dg  or  alternatively,  a 

thick  parallel  wake;  and  (c)  flow  separation  from  the  base  and  re-attachment  to  a thin  sting  of  diameter  dg  / 2. 

These  examples  give  a wide  variation  in  base  pressure,  but  the  influence  of  the  base  flow  extends  only  about  one 
base  diameter  upstream  from  the  base,  giving  confidence  in  the  validity  of  the  solution  ahead  of  this  region. 

If  condition  (c)  of  figure  2 is  a fair  representation  of  the  base  flow  with  a thin  sting  present,  then  condition  (b) 
implies  that  a large  sting  can  cause  a significant  strengthening  of  the  adverse  pressure  gradient  near  the  model  base 
and  perhaps  influence  the  measured  fin  characteristics.  The  uncertainty  regarding  the  best  choice  of  base  flow 
geometry  can  be  partly  resolved  from  measured  base-pressure  data.  Such  data  are  available  from  reference  S for 
a variant  Afi)  of  the  Body-A  configuration,  which  is  discussed  in  the  following  section. 
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3.2  Conriguration  A(i) 

The  A(i)  variant  has  the  rear  portion  of  the  standard  boattaii,  where  fins  are  normally  located,  replaced  by  a 
cylindrical  section  which  gives  a 45%  increase  in  d_  over  the  standard  conftguration,  the  slope  of  the  tail  cone 

D 

section  remaining  the  same  at  7°  40''  . The  afterbody  of  this  configuration  and  the  calculated  pressure 
distribution  are  shown  in  figure  3(a),  with  the  flow  assumed  to  separate  from  the  base  parallel  to  the  cylindrical 
section.  The  computed  values  not  shown  for  the  forebody  are  identical  with  those  given  in  figure  2 for  the 
standard  configuration.  Figure  3(a)  shows  a strong,  favourable  pressure  gradient  along  the  tail  cylinder,  instead 
of  the  weak  adverse  gradient  on  the  standard  body  (figure  2 condition  (c)).  This  may  result  in  improved  fin 
lifting  performance  by  reducing  the  tendency  for  flow  separation  in  the  lee-side  fin  root  region. 

Figure  3(b)  shows  the  effect  of  removing  the  surface  slope  discontinuity  on  the  A(i)  afterbody  by  fitting  a cubic 
polynominal  to  the  conical  and  cylindrical  sections  between  X = lOandx  = 11.  Similar  fairing  will  be  effected 
by  the  boundary  layer  in  the  real  flow,  and  wind-tunnel  measurements  (from  reference  5)  on  the  A(i)  model, 
which  are  included  on  figure  3(b),  compare  favourably  with  the  present  results.  A sting  of  diameter  ^ 

for  the  measurements,  and  the  data  indicate  some  flow  expansion  behind  the  base.  The  calculated  results  can  be 
made  to  agree  more  closely  by  including  a small  expansion  in  the  base-flow  model,  as  shown  in  figure  3(c). 

3.3  Configuration  A(ii) 

The  A(ii)  configuration  has  a flared  tail  section  instead  of  the  cylindrical  section  of  the  A(i)  variant.  This  is 
shown  in  figure  4 together  with  the  strong,  favourable  pressure  gradient  resulting  from  this  geometry.  As 
mentioned  above  this  is  expected  to  be  beneficial  to  fin  performance,  but  to  offset  this  there  is  a strong  adverse 
gradient  ahead  of  the  flared  section  which  is  likely  to  cause  some  degree  of  flow  separation,  and  higher  drag 
resulting  from  the  larger  base  area. 

3.4  Body-A  with  variable  afterbody  slope 

Pressure  distributions  were  calculated  on  the  Body-A  forebody  fitted  with  conical  afterbodies  of  half  angle 
5°,  10°,  15°  and  20°  respectively.  The  results  are  displayed  in  figure  5.  The  adverse  pressure  gradient  on  the 
afterbody  strengthens  dramatically  with  increasing  slope,  but  the  afterbody  slope  has  little  effect  on  pressure 
coefficients  ahead  of  about  the  mid-point  of  the  central  cylindrical  section.  If  it  were  desired  to  find  the  maximum 
slope  tolerated  before  separation  became  significant,  a boundary-layer  program  could  be  used  in  conjunction  with 
the  present  analysis. 

3.5  Body-B  results 

The  relatively  small  base  of  the  Body-B  configuration,  shown  in  figure  6, would  require  the  sting  diameter  for 
small  wind-tunnel  models  to  be  as  great  or  greater  than  d_  . Figure  6 compares  the  pressure  distributions  near 

D 

the  base  when  a sting  of  diameter  d„  is  present  and  absent.  Consistent  with  the  earlier  results,  the  sting  magni- 

D 

fies  the  unfavourable  pressure  gradient  over  the  tail  section  where  fins  would  normally  be  located.  This  should 
be  considered  when  interpreting  measurements  of  fin  properties. 


4.  CONCLUSIONS 

Pressure  distributions  have  been  calculated  on  two  bodies  of  revolution  in  incompressible,  inviscid,  axisymmetric  flow 
The  programs  of  Albone(ref.2)  give  sufficiently  accurate  results  with  minimal  user  effort  and  computing  time.  Bodies 
with  bluff  bases  such  as  boattails  can  be  treated  by  assuming  the  body  surface  to  be  extended,  thereby  simulating  the 
flow  separation  from  the  base.  The  choice  of  an  appropriate  base-flow  model  can  be  made  more  reliable  if  base- 
pressure  measurements  are  available  for  comparison. 

The  results  have  shown  that  a favourable  pressure  gradient  can  be  created  by  replacing  the  rear  portion  of  a boattaii 
with  a cylindrical  section,  which  may  improve  the  effectiveness  of  fins  mounted  thereon.  Conversely,  fin  lift  measure- 
ments on  wind-tunnel  models  may  be  adversely  effected  by  the  presence  of  a sting  if  its  diameter  is  comparable  to  the 
model  base  diameter. 
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(a)  Closed-body  program 
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C PHOC'W.M'  FflR  Cl.nSFn  AX  ISYr'MFTWIC  RnOY 

n lA'FNF  I ON  A(A0),X<AP),F(4P),Fl(40)tF(40).Fl(40).F?>(40).Xl(A0).U(A? 
I ) .CP(4?  ) ,XK  (40.40)  .F2(40) 

C N'MAX  IS  MAX  MO  OF  ITFRATIOMS  ALLOWFO.  FPS  IS  ARSOLUTF  FPPOR  LIMIT 

C BFTWFFM  iterations  at  CONVFROFNCF.  XXO.XXl  ARE  FORF  AMD  AFT 

C AM5CISSAF  OF  fAODY.  S/R  f>Y  OFMFRATFS  ORDINATE  FB  AND  SLOPE  FIP  OF 

C BODY  AT  EACH  ABSCISSA  XB . 

REAL  MACH 

DATA  NMAX.FTJS/  SO. 0.0001/ 

DATA  XXO . XX 1 /O .0 . I .0/ 

DATA  .VACH/0.0/ 

RFAL»H  ADUNMPC)  /.004S?1?771 ..0104R82B4S..01642105fl4..022P4SP49?. 

1 .0?7o.T70070..0334S01or'S..07B7B21SP..04:iB7090B?.  .04P6RSP0  76  . 

2 .053227P47..05743P7pai..061306242S..Ofi4P040135..0S79120458. 

3 .0  70B 1 1 64  74  . .072PP6SP24 .,0747231691  ..0761103619.  .0770398182 . 

4 .077S0S948C/ 

RF  ALPPyniiv  ( 20  ) /. 99823771  ..900726239  . , 9772S99S  . . 9S79  1 68  1 9 . 

1 . 932  M 2 PC'H  ..  90  2 0 9880  7 . .M6S959S0  3 . .824  61223  1 . . 778  30S6S  1 . 

2 . 7273  1 82 '■.B,.6719S66f'S..612SS3H9,. 549467125. .4  83078802. .413779204 

3 . .34  1 >8-/)C'91  . .2681  52  1 P5  ..  1 9269758  1 . .1  1 60  R4  07  1 . .038772418/ 

Bf-  TSO-  1 -rACn«yi>CH 

BF  T A r.SoiJ  T ( Ml  - TSD  ) 

XM 1 0=  (>;x;  +XX 1 ) /2  . 

CALL  Kf'DY  ( X>' ID  . YK  1 0 ,F  1 B ) 

YM  I D = Y'-  I D>:<BFTA 
F I B=F  1 M<<nFTA 
DO  2 K=I .20 

KK  =4 1 -X 
A ( K 1 ^Apot-'  ( K ) 

A (KK  ) = A (K  1 

X ( KK  ) =Xni'M  ( K ) 

2 X (K ) =-xnHM(K 1 
no  8 K= 1 ,40 

X(K1=C).S4(X(K14(XX1-XX01  + XaitaX0) 

XB=X (K 1 

CALL  BPOY ( XR ,FP ,F 1 B ) 

Fn=FR  *BF  TA 
F I B=F  1 B-EI'FTA 
F (K I =FB 
F 1 ( K ) =F  1 B 
F2 (K 1 =F IB 

XI  (K  ) = (X(K)-XX0)1'(XX1-X(K)  1 

FI  (K  1 = 1 ./SORT ( 1 . +F 1 (K  ) 4F 1 ( K 1 1 

8 F <K  ) =F (K  ) *F ( K ) 

XK0=XKI ( <XX1-XX0>/(2.*YM1D) ) 
on  3 1=1,40 

FOFT=F ( I )/X  I ( I 1 
XL  = 50RT ( 1 ./FOFT 1 
XK2=(1.+XK0)/{1.+XK1(XL)1 
SUM  = 0 . 

no  4 j= 1 ,4 0 

XDIF=(X(J)-X(I))4(X(J)-X(I)) 

XK  ( 1 , J 1=F ( J 1 /SORT ( (XDIF  + F(J)  )A43l 
r,XT  = FOFT«Xl  ( J ) 
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XKnASn-c',XT/snWT((XDIF  + GXT)’>*3) 

A Pl)M=Snr.’  + A(J)«(XK(  I ,J)-XKDASH) 

F(  I )=1  1 . + XK0)’MXX1-XX0)  FSUM-XK? 

F1(I)=(F(I)+XKC)/(1+XK0) 

3 IK  I ) = ( 1 +XK0)«F1  ( 1 ) 

N = 0 

FMAX=1 OC 
r-,  N = N+l 

FLAST=FMAX 

nn  fi  I =1 , AO 
U ( I ) =U  ( 1 ) +F  1 ( I ) *E  ( I > 

SUM=0 . 

nn  7 j=i,Ao 

7 SIIM=SUM  + A ( J ) «XK  ( 1 , J ) * ( E ( J ) -F  ( I ) ) 

F>  F2(I)=F1(I)*F(I  )-  .25*SUM«  ( XX  1 - XXO  ) 
on  13  K = 1 , AO 
13  F (K  ) =F2 (K  ) 

FMAX  = AHS (F { 1 ) ) 

0012  K=1 ,A0 

IF(AnS(F(K)).r,T.FMAX)EMAX  = ARS(F(K)) 

1?  CONTINUE 

IF(FLAST-FMAX)  lA,17tl7 

17  IF (fkax.gt.fps.and.n.ut.nmax ) on  TO  5 

IF  WU I TF  ( 6 , 1 A )N .FMAX  .MACH 

lA  Fm?MAT(20H  NO.  OF  ITERATIONS  =.I3»10X.13H  MAX.  ERROR  =.F15.7. 

1 lOX.QHMACH  Nn.=.F5.3//) 
no  P K=1 . AO 
F (K  )=SnRT(F  (K  ) ) 

F ( K ) =F (K  ) /HF  T A 
F2 (K  ) =F2  <K  ) /npT A 
IKK  >=1  +(U(K  )-l  l/OETSO 
IF(MACH-O.l)  2R.2R.1R 
2<-i  C IMC  = I -IK  K ) >KK  K ) 

CP(K)=CINC  + 0.?S>1<MACH>I'MACH'«CINC>!‘CINC 
r.n  TO  “ 

IQ  CP(K  ) = ( ( 1-0.2’>MACH*MACHA(IKK)>KU(K)-1  ) )X>*3.5-1  )/(0.7*maCH*MACH) 

P CONTINUE 
WW I TF ( A , 1 5 > 

IF  FORMAT(AX.?H  X.13X.2H  Y.12X.AH  OY/DX.TX.RH  VE  LOC  I.  TY  . 3X  . 1 5H  PRESS. 
ICOFFFT.  ) 

WRITF(A,1A)(X(K),F(K).F2(K).U(K),CP(K).K=1.A0) 
ir,  FPPI''AT(3(F10.A.AX).Ffl,A.6X.FR.A) 

■■<7  STOP 
FMO 

FUMC  T I ON  XK 1 ( n ) 

IF  ( M- 1 .000 1 ) 3 . 3 . 1 
7 IF  ( F'-O  . POPP  ) 2 . A . A 
A XK  1 =0  .o 
RF  TURN 
1 C=0*O 

n = SOPT (C-1  . ) 

E = ALOfi  ( M +n  ) *n 

XK 1 = ( F-n ) / ( C«0-F  ) 

nr  TURN 

P C-M70 

o-soPT ( 1 .-r  ) 

F =Al.or,  ( ( 1 ,+n  ) /FI ) *C 

XK  i = (o-F  ) / (?'  ..i-n^oAO-D+E  ) 

Ri-  TURN 
FMO 


non 
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(b)  InPinite  ■ afterbody  program 


C PROGRAM  FOR  BOOV  FXTFNOIN'G  TO  INFINITY 

O  IMF  NS  ION  At  F.0  ) . X ( fiO  ) ,F1  ( ftO  ) .F  ( fiO  ) .F  1 < 60  ) ,F?  ( fiO  ) .U(  60  ) .F  (60  ) tCP(  ftO 
1 ) .F? (60  ) .XK  (60  t60 ) 

NMAX  IS  MAX  NO  OF  ITFRATIONS  ALLOWFO.  CONVFRGFNCF  WHFN  MAX  DIFFFRFNCF 
BFTWFFN  ITFPATIONS  IS  L.F.  FPS.  XXO  IS  ABSCISSA  OF  NOSFtN.B.  JUNCTION 
OF  BOOY/AF TFRBODY  IS  AT  X=0).  YAFT  IS  RADIUS  OF  PARALLFL  AFTFRBODY. 

C S-R  nnOY  <,FMFRATFS  radius  fb  and  SLOPF  fib  at  FACH  abscissa  xb 
REAL  I.FNGTH 
RFAL  MACH 

DATA  NMAX, FPS/  50,0.0001/ 

DATA  LFMGTH.YAFT  /2.26B,0.5/ 

DATA  MACH/0.0/ 

RR  AL’l'H  AOUM  (15)  /.007Q6BlO?5,.01B4  664  6B3,  ,02B784  707R  . .0387991926  . 

1 . 04 H40267PB , .0574931 562 , .0659742299, .0737559747 , .0 BO 7558952 , 

2 . 0B6BOO7B7? ,.0921225222, .0063687372  , .0995934206, .101762390, 

3 .102^52653/ 

RF  AL««  XOUM  ( 15  ) / .906  80,34  84  , .08  3668  123,  .06002  1865  , . 92620004  7 , 

1 .HF2560S36, .820565762 ,.767777432,. 607850495 ,.620526183, .536624148 

2 . .447  0 3377 ,.  35  2 704  726  ,. 254636926 , . 15  38690  14 ,,05147184  26/ 
BFTSO=l-f'ACH<<MACH 
MrTA  = SORTCiFTSO) 

YAFT  = YAFT‘^nFTA/LFNGTH 

XXli=-l  , 

on  2 K ^1  , 1 5 
KK=31 -K 
A (X ) =AOU” (X ) 

A ( X X ) = A ( X ) 

X ( X 1 =xni  iM  ( X ) 

X ( X K ) - X ( X ) 

2 X (K  ) --X  (X  ) 
nn  17  K = 1 , TO 
XX  =30  +K 
A ( XK  ) =A  ( X ) 

17  X (KK  )=X (K  ) 

nn  8 K=1 ,30 

X(K)=0.5*XX0*(1.-X(X)) 

XHt  lfng  TH<=  ( X ( K ) + 1 . ) 

CALL  BODY ( XB ,FB,F 1 B ) 

fb=fb/  length 

FB=FB«Rft A 
P1R=F  lM*riFTA 
F (K ) =FB«FB 

F1(K)  = 1./  SORT(  I .-FFIB-XFIP) 

F2(K )=F1B 

8 A (X )=0 .25«XX0*A (K 1 

nni4K=31 ,60 

X(K)=(1.+X(K))/(1,-X(K) ) 

F(K )=YAFT«YAFT 
FI (K )=1 , 

F2(K  )=0  . 

14  A(K)=-0,25*(1.+X(K))4(1.+X(K))4A(K) 

nn  3 1=1,60 
F0FT=F( I)/(X(I)-XX0) 

F(  I )=0. 

^ ^ M .»  " ■ ^ 

f '-i  V * Wi  '•  - 


nr.c 
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no  4 j = i ,fso 

Xr>IF={X(J)-X(I  ) )4(X(J)-X(I  )) 

XK(r.J)=F(J)  / SORT(  (XOIF  + F(J  ) >443) 

CXT=FPFT4(X(J)-XX0) 

XKPASH  = GXT/(  XntF+r,XT  >441 .5 

4 F ( 1 ) =F ( I ) + A ( J > 4 ( XK ( I , J )-XKOASH  ) 

FI  ( 1 )=F ( I ) 

1 U(  1 )=F1  ( I ) 

N = 0 

5 M=M+I 
EMAX=0 
DOft  1=1,60 

l)(  1 ) =U  ( I)  +F  1 ( 1 > 4F  ( 1 ) 

F?(  1 )=FI  < I )4F(  1 ) 

DO  7 J=1 , 60 

7 F?(I)=F2(I)+A(J)4XK(l,J)4(e(J)-F<I)> 

6 CnN  T IMUF 

no  IP  1=1,60 
F ( 1 ) =F? ( I ) 

1?  IF(  AOS ( F ( 1 ) ) ,6T  ,FMAX  ) FMAX=  ARS(F(I>) 

IF  ( FM  AX  ,GT  ,FPS,  AMO  ,N  ,LT,NV<AX  ) nn  TO  S 
'VI' 1 TF  ( r.  , 10  )N  ,FM  AX  .MACH 

10  rr!OMAT(?0H  ^'n,  OF  ITFRATIONF  =,13,10X,13H  MAX.  FPROR  =»F16.7, 

1 lOX,  IOHMaCH  M0.=  ,F4.?,//) 

Ul  IM  I 1)  K -•  1 , 6(1 

I-  (K  soft  (f  (K  ) ) 

F r K ) =LF  N6TH4F  < K,  ) /OF  TA 
X ( K ) =LrM6TH* ( X ( K ) + 1 . ) 

FP(K  )=F?(K  l/OFTA 
IKK  ) = 1 + (I)(K  )-l  )/HFTSO 
IF  (IWACH-0 . 1 ) ?o,?R.lR 
/o  r.  I MC  = 1 -II  (K  ) *U(K  ) 

Cl>  ( K ) =r  IMC +0  . ?6  4M  ACM  4MACH4C  I MC4C  I NC 
60  TO  O 

1 O ri->(K  ) = ( ( 1-0.?4MACH4MACH4(u(K  )4U(K  1-1  ) >443.6-1  >/(0.74MACH4MACH> 
o CONTI  Ml  IF 

■VRI  TF  ( 6 , 1 6 > 

16  FOR"AT(6X,?H  X,13X,?H  Y,1?X.6H  DY/DX.7X.QH  VFLOCITY.RX, 

1 14H  PPFSS.  COFFT. ) 

'"WlTr(6,16>(X(K),F(K>,F2(K>,O(K>,CP(K>,K  = l,40> 

16  FnRMAT(F10.6,6X,F10.6,6X,F10.6,6X,Fn.4,6X,F«.4> 

r 

/i-^  STOP 
FNO 


DccT  rnr/ 

^ w t > 
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Figure  3(a),  (b)  and  (c) 


(a)  Basic  cone-cylinder 


(b)  With  fairing 


30  40  60  60 

x(CKl) 

(c)  Modified  base  flow 


Figure  3.  Pressure  distributions  on  A(i)  afterbody 


Figure  4.  Pressure  distribution  on  Body-A(ii) 


(b)  Sting 
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